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ABSTRACT 

The Array for Microwave Background Anisotropy (AMiBA) is a radio interferometer for research in cosmology, 
currently operating 7 0.6m diameter antennas co-mounted on a 6m diameter platform driven by a hexapod 
mount. AMiBA is currently the largest hexapod telescope. We briefly summarize the hexapod operation with 
the current pointing error model. We then focus on the upcoming 13-element expansion with its potential 
difficulties and solutions. Photogrammetry measurements of the platform reveal deformations at a level which 
can affect the optical pointing and the receiver radio phase. In order to prepare for the 13-element upgrade, two 
optical telescopes are installed on the platform to correlate optical pointing tests. Being mounted on different 
locations, the residuals of the two sets of pointing errors show a characteristic phase and amplitude difference 
as a function of the platform deformation pattern. These results depend on the telescope's azimuth, elevation 
and polarization position. An analytical model for the deformation is derived in order to separate the local 
deformation induced error from the real hexapod pointing error. Similarly, we demonstrate that the deformation 
induced radio phase error can be reliably modeled and calibrated, which allows us to recover the ideal synthesized 
beam in amplitude and shape of up to 90% or more. The resulting array efficiency and its limits are discussed 
based on the derived errors. 

1. INTRODUCTION 

The Array for Microwave Background Anisotropy (AMiBA) is a forefront radio interferometer for research in 
cosmology. The project is led by the Acadcmia Sinica, Institute of Astronomy and Astrophysics (ASIAA), 
Taiwan, with major collaborations with National Taiwan University, Physics Department (NTUP), Electrical 
Engineering Department (NTUEE), and the Australian Telescope National Facility (ATNF). Contributions also 
came from the Carnegie Mellon University and NRAO. As a dual-channel 86-102 GHz interferometer array of 
up to 19 elements, AMiBA is designed to have full polarization capabilities, sampling structures greater than 
2' in size. The AMiBA targets the distribution of clusters of galaxies via the Sunyaev-Zel'dovich (SZ) Effect 1 . 
It will also measure the Cosmic Microwave Background (CMB) temperature anisotropics 2 on scales which are 
sensitive to structure formation scenarios of the Universe. AMiBA is sited on Mauna Loa, Big Island, Hawaii, 
at an elevation of 3,425m, currently operating the initial phase with 7 antennas of 0.6m diameter in a compact 
configuration, FigfT] The correlator and the receivers with the antennas are co-mounted on a fully steerable 6m 
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Figure 1. Front view of the AMiBA. Installed are 7 antennas in compact configuration, giving 0.6m, 1.04m and 1.2m 
baselines. Free receiver holes in the platform allow for different array configurations and for the expansion phase with 
13 antennas. The 1st 8 inch refractor (OT1) for optical pointing is attached to the black bracket below the lowermost 
antenna at a distance of 1.4m from the platform center. A second optical telescope (OT2) is located at 90° to the left 
hand side at the same radius. 



diameter platform controlled by a hexapod mount. The AMiBA hexapod mount with its local control system 
was designed and manufactured by Vertex Antennentechnik GmbH, Duisburg, Germany. The 0.6m diameter 
Cassegrain antenna 3 field-of-view (fov) is ~ 22' at Full Width Half Maximum (FWHM). The synthesized beam 
resolution in the hexagonally compact configuration (0.6m, 1.04m and 1.2m baselines) is ~ 6'. This sets our 
target pointing accuracy to about ~ 0.6'. 

In parallel to the current 7-element observation program, AMiBA is undergoing a 13-element expansion with 
1.2m diameter Cassegrain antennas (fov ~ 11'). The resulting synthesized beam of about ~ 2' requires then 
ideally an improved pointing accuracy of about ~ 0.2'. The platform photogrammetry tests 4 ' 5 have revealed 
local deformations which might affect individual antenna pointings (including the optical telescope pointing) and 
their radio phases for the longer baselines. It therefore seems desirable to include such a correction in the current 
pointing error model and in the phase correction scheme. 

The goal of this paper is to investigate a refined pointing error and phase correction model for the AMiBA 
expansion phase. The adopted approach can be of general interest since a platform deformation for a hexapod 
of this size is likely to be a generic problem depending on the level of accuracy which is needed. The paper is 
organized as following: Section [2] summarizes the pointing error model used for the currently operating 7-element 
compact configuration. Section [3] presents the approach taken to improve the pointing accuracy and the phase 
correction for the 13-element expansion and it also discusses its limitations. Previous AMiBA status reports and 
the current observations are described in a series of papers. 6,7,8 ' 9 

2. 7-ELEMENT POINTING 

The AMiBA hexapod has an elevation limit of 30° with a continuous azimuth range and a polarization limit 
of ±30° at each azimuth-elevation position. An important consequence of the hexapod mount are the non- 
existent azimuth and elevation encoders compared to more traditional telescopes. The 3-dimensional locations 
of the upper and lower universal joints in the reference positions and the 6 exact encoder-measured jack lengths 
completely define the geometry of the mount. The 6 upper universal joint locations define a best-fit plane with 
its normal defining the resulting pointing axis. During the operation it is thus most essential to control the 
lengths of the 6 jacks. In order to achieve a pointing accuracy of ~ 0.6', a 2-step approach 4 was adopted: 



Figure 2. OT1 (in the back) and OT2 (in the front, lower left hand side corner), mounted with a 90° phase shift in between 
them at a radius of 1.4m. Clearly seen are also the reflecting photogrammetry targets. 



• first pointing run: 

all known pointing corrections are activated. The star position on the optical telescope (OT) ccd is split 
into an azimuth and an elevation error as a function of the mount position. The errors are analyzed to 
separate the OT collimation error from the mount pointing error. The resulting residual errors are filled 
in a 3-dimensional error interpolation table (IT) as a function of elevation, azimuth and polarization. 

• second pointing run: 

all known pointing corrections and the IT are activated and the improvement is verified. 

The pointing error model can be listed in two groups: 

• Jack corrections: 

- jack pitch error compensation: 

calibrated at reference temperature, look-up table as a function of jack length for each jack 

- jack temperature compensation: 

analytical model, real-time measurement with 3 temperature sensors along the jacks 

- jack rotation correction: 

analytical model for encoder undetected jack rotation 

- support cone compensation mode: 

model for universal joint shifts due to temperature changes, monitored with several temperature sensors 

• Telescope corrections: 

- radio/optical refraction: 

analytical refraction model, wavelength dependent 

- OT correction (collimation error): 

analytical model, collimation correction derived from pointing data 

- anchor cone orientation correction: 

analytical model, rotation of mount with respect to north 

- error interpolation table (IT): 
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Figure 3. Left hand side figure: Polar plot of the total pointing error (azimuth and elevation component combined) without 
IT correction. The OT collimation error is shown as a circle which represents the OT tilt ~ 9' with respect to the mount 
pointing axis. The data points scattered off the circle show the residual pointing errors for the IT, leading to ~ 1' rms 
pointing error. The test was done with 250 stars. Right hand side figure: Polar plot of the total pointing error with the 
IT correction. The reduced scatter verifies the improved pointing with ~ 0.4' rms pointing error. 



interpolation table, residual errors which are not explicitly modeled, derived from pointing data 



The 2-step approach with the above outlined pointing error model and the additional improvement from the IT 
is illustrated in FigJ3] The original pointing error from ~ 1' rms is reduced to ~ 0.4' rms with the help of the 
IT. These first results were derived with a single OT (OT1) on the platform (Fig(2). 



3. 13-ELEMENT POINTING 

3.1 Comparison OT1 - OT2 

The mount pointing can only account for an average pointing. Any possible local irregularity or (platform) 
deformation is averaged out and therefore not taken properly into account. This is in particular also the case for 
each antenna and the OTs. Therefore, the results in section [5] have to be re-analyzed for an improved pointing. 
A second optical telescope (OT2) was installed on the platform in order to cross-check the pointing results from 
OT1. This also provides an independent and complementary tool to the photogrammetry tests. The OT2 (SBIG 
ST-237A ccd with a hyperstar f/2.0. lens, fov ~ 42' x 32') is installed at an angle of 90° with respect to the 
OT1 at the same radius r=1.4m from the platform center, Fig[2] It is also integrated in an automatic pointing 
schedule 4 which allows us to simultaneously take ccd images with both OTs. 

The Figsf4] and [5] show the separated azimuth and elevation errors for both OT1 and OT2 after the OTs 
collimation errors (solid circle in Fig[3]) have been removed. The pointing run was done with 200 stars, equally 
distributed with constant solid angle over the total accessible sky. The platform polarization was kept at zero for 
this test. The errors are then interpolated (cubic interpolation) on a regular azimuth-elevation grid. Apparent 
similar pointing errors are measured by both OT1 and OT2, Figs|4] and E] upper two panels, for both azimuth 
and elevation errors with some clear dominating patches at lower elevation. The overall azimuth and elevation 
rms errors between OT1 and OT2 agree within 4 — 7". 

However, taking the difference between OT1 and OT2 for their corresponding azimuth and elevation errors 
reveals clear systematic patterns (Figs0] and [3] lower two panels) with an amplitude of up to 1'. 

In the ideal case, this difference should be negligible. If there is a local error contribution at each OT location, 
this measures the difference between the local errors, symbolically written as: 



error ot\ — localori + pntgerror 
error QT2 = localoT2 + pntgerror 
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Figure 4. Upper two panels: OT1 and OT2 azimuth errors after removing the OTs collimation errors. Both OTs show 
similar dominating structures. Lower two panels: difference in azimuth errors shown in a contour plot and extracted for 
some elevation bands. A clear oscillating structure with positive and negative domains becomes visible. All units are in 
arcmin. 
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Figure 5. As in Fig(4]but for the elevation error. Again, a similar structure in the OT1 and OT2 data is apparent, as well 
as the oscillating feature. 
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Figure 6. Left hand side figure: A typical saddle-type platform deformation in z-direction, measured from photogrammetry 
data. The units are in mm. The two black dots along the black circle mark the locations of the two OTs. Right hand 
side figure: Platform deformation in z-direction as a function of platform azimuth at a radius r=1.4m where the OTs 
are mounted. The selected data points are within an annulus of r=1.4m±0.2m, extracted along the black circle in 
the left hand side figure. The maximum deformation is measured at the extreme position (az,el,pol) = (0,30,20) and sets 
therefore an upper limit to the deformation error of about 1'. Also illustrated for comparison is an average mount position 
(az,el,pol) = (0,60,00) with an amplitude of about 0.11 mm, leading to a corresponding OT tilt of about 0.25'. Clearly 
visible is the saddle structure with a functional form cos(2az). 



dif ference = localon — localoT2 



The difference plots show clear 'plus' and 'minus' patches, which are likely the result of the platform deformation 
(section !3-2[) which sometimes tilts the OTs towards each other and sometimes away from each other. Over a 
full azimuth circle the OTs are aligned with each other four times (2 local maxima and 2 local minima on the 
saddle pattern, FigJS]) and the difference becomes zero. The bottom panels in the FigsH] and [5] illustrate this for 
some elevation bands. 

3.2 Platform deformation modeling 
3.2.1 Optical pointing 

Since the mount pointing can only correct for an average pointing error, any local deformation error should 
be separated. In order to do that, the photogrammetry results 4 are used to construct an analytical model for 
the local platform deformation induced pointing error. The saddle-type deformation in z-direction at a fixed 
telescope position can be modeled fairly well with a function of type f(x,y) = A(x 2 — y 2 ), both for the entire 
platform and for an annulus at a radius r (Fig[6]). The OT at any location has its optical pointing axes normal 
to the tangential plane at the deforming platform. The function / describing the deformation with its gradient 
V/(r, a) — in spherical coordinates x = rcosa, y = rsina are written as: 

Ar 2 (cos 2 a - sin 2 a) , ( 1 ) 

2Ar(cos 2 a-sin 2 a), (2) 

— 4Arcosasina. (3) 

In order to calculate the azimuth and elevation error components for the induced pointing error, the angles 
d between the normal of the tangential plane (h) and the normal to the platform no along the azimuthal and 
radial direction are calculated as a function of the platform azimuth a and radius r: 



f(r, a) 
df 
dr 
10/ 
r da 



d = Z(n ,n) = 7r/2 - Z(n"b,t), 



(4) 



with t being either one of the unit tangential vectors derived in Eqs.([2]) and ©. Using the scalar product 
cosG?(no, fi) — n@ • h, we find the deformation induced pointing errors in azimuthal and radial direction as: 



where r = R is the radius where the OTs are mounted. Having two OTs mounted on the platform at the 
same radius R with a separation of 90° in between them (Fig[2|) provides for each of the error components in 
Eqs.f5|) and (J6j) a set of two equations with a characteristic 90° phase shift in a. Additionally, the saddle pattern 
has an offset (j> relating the saddle pattern orientation to the mount position, i.e. a — az + <f> for OT1 and 
a = az + 4> + it/2 for OT2. The deformation is then entirely described with the functional form in the Eqs.([5]) 
and ([6]) and the two fitting parameters A and <f>. 

We remark that the Eqs.([5]) and ^ are almost identical to the function of type Acos(2az + <f>). For simplicity 
we use the later one in section [5751 where A and R are absorbed in the amplitude A. 

3.2.2 Radio pointing: radio phase correction and synthesized beam 

Besides the OTs, the platform deformation also affects each individual radio antenna: the tilting of each an- 
tenna leads to a changing misalignment error and therefore reduces the efficiency (section 13. 4p . the deformation 
correction in z-direction gives a radio phase error. Like the OTs measure directly the relevant quantity for the 
pointing error analysis (error in arcmin), the photogrammetry directly gives the relative phase errors (in mm) 
with respect to a reference position for all antenna locations. Fig[7] shows the distorted synthesized beams for a 
set of mount positions with their measured radio phase errors. For the largest deformations (up to 800/zm rms), 
the beam for a 13-element configuration is severely distorted. Since the platform deformations are repeatable 
(within ~ 50/im), the deformation can be modeled. Based on a set of about 50 measured mount positions, we 
have developed both an analytical modeQ and a 3-dimensional interpolation table approach in order to predict 
and correct the phase error. The deformation rms residuals after correction are less then 100/im, which restores 
the synthesized beam amplitude to more than 90%. For the current 7-element compact configuration the beam 
distortion is negligible. 



3.3 Correlation OT1 - OT2: local platform deformation and remaining pointing error 

The goal in this section is to separate the real mount pointing error from the local platform deformation which 
affects the two OTs. This will then define the pointing accuracy. In order to achieve this goal, we combine the 
data from OT1 and OT2, each of them providing an azimuth and an elevation error for each telescope position. 
We assume that both OTs are measuring the same underlying pointing error (p az in azimuth and p e i elevation) 
which is perturbed by the local deformation (d az , d e i) which we describe with the analytical model derived in 
section 13.21 A simultaneous fitting (non-linear curve fitting optimization in the least-square sense) is done for 4 
unknowns (deformation amplitude A, deformation phase <fi, azimuth pointing error p az , elevation pointing error 



d, 




(•5) 



* The model is again based on the functional form Acos(2az + <f>), with the amplitude A and the phase <f> being 
complicated functions of azimuth, elevation and polarization. 
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Figure 7. Top panels: The distorted synthesized beams (labeled as 'Image') for the 7-element and 13-element configuration 
with their best-fit Gaussian beams. Different curves are labeled with the mount position in azimuth-elevation-polarization, 
with the first curve being a reference position without deformation. An 'Image' is calculated from the measured phase 
errors and represents therefore a point source image. For the current 7-element compact configuration (left hand side 
panel) the distortion is negligible. Bottom panels: Gaussian beam peak and FWHM ratios (with respect to an ideal 
no-deformation case) for different configurations at different mount positions. The ratios decrease with lower elevation. 
1.4+2.4m/1.4-|-2.8m are different possible configurations for 13-element. 




Figure 8. Left hand side: Original platform deformation rms and residual rms after model correction for a data set of 
about 40 telescope positions. Right hand side: S re i and 5 a bs (in long-dashed lines) for the efficiency calculations. The two 
antenna primary beams and the synthesized beam are shown in solid and short-dashed circles, respectively. 



Pel ) with a system of 4 coupled equation^] 



Paz 1 = d az {az,A,4>) +p az , (7) 

P°J 2 = d az (az,A,<l>+ ~)+p az) (8) 

P ^ T1 = d el {az,A,<j>)+ Peh (9) 

P°J 2 = d el (az,A,<f>+ |)+p d , (10) 

where Pa-f 1 , P°I 2 > P^ 1 an d P^ T2 are the errors measured on the CCDs after the OTs collimation errors are 
removed (FigsH] and [5]). az is the mount azimuth position for each set of pointing data, d is the function 
describing the platform deformation, derived in section 13.21 The deformation terms d for OTl and OT2 are 
coupled through a phase shift 7r/2, because the two OTs are mounted on the platform with a 90° separation in 
between thenljj . 



Fig[5]shows the result for the 4 fitting parameters. Even after removing the platform deformation component, 
the azimuth and elevation errors show similar dominating structures as in the original data sets from OTl and 
OT2. The overall positive azimuth error is likely due to a remaining cone mis-orientation with respect to north. 
The existence of this real dominating structure also explains why the deformation corrected pointing errors are 
not random. The newly derived values for p az and p e i are then used for the IT. The typical local deformation 
introduces azimuth and elevation pointing errors between —V and +1'. The deformation corrected pointing 
errors are between —1' and up to 2 — 3'. 

As a consistency check we show the separated azimuth and elevation deformation errors for OTl and OT2 in 
FigHOl The difference between the local deformation errors between OTl and OT2 are then to be compared 
with the difference between the raw OTl and OT2 errors without any modeling (error separation) in the bottom 
panels in the Figs 0] and E] The difference plots indeed agree very well. 

3.4 Pointing limit, phase error and array efficiency 

The platform deformation (and any other local error) complicates the definition and the extraction of the mount 
pointing error. The question rises how far the pointing accuracy can be improved before the local deformations 
start to dominate. In order to answer this question we quantify the array efficiency loss resulting from the mount 
pointing error and the local platform deformation. For completeness we also include a radio misalignment error 
for each baseline (each antenna pair). For each antenna pair we calculate a relative and an absolute error, 5 re i 
and S ci bs, respectively. S re i measures the shift in the overlap of two antenna primary beams (relative antenna 
misalignment) due to any local irregularities, whereas 5 aos quantifies the (absolute) synthesized beam shift for 
each antenna pair with respect to a perfect pointing axis (Fig[8l right hand side): 



Srel = \J {{d\MZ - d 2 .az) + (pV^az - m 2 ^z)) 2 + {{d\,el - d 2 ,el) + (™l,ei ~ ™2,ei)) 2 > (H) 

, fdi a z + <k,az , mi,az + m 2 ,a,z , \ 2 , { di e l + d 2 , e l , mi,el + ni 2 ei \ 2 
Sabs = \ = + = + Paz) + = + = + Pel) > (12) 



where dk and nik are the platform deformation and the misalignment error, respectively, for k = 1,2 for each 
antenna pair, each having an azimuth and elevation error component. p az and p e i are the mount azimuth and 



' Separating platform and pointing error at each mount position (with only 4 measurements for 4 equations) has been 
found to be more reliable than trying to extract a deformation correction for an entire mount azimuth or elevation range. 
This is due to the fact that for a single specific mount position the deformation at the OTs locations can very well be 
modeled by a function of type A cos(2az + (j>) (Fig©, whereas for an azimuth or elevation range both phase and amplitude 
are complicated functions showing large variations (reflected in the bottom panel in the Figs |4] and [5|. 

* For practical purposes, we use d az = j4cos(2az + 0) and d e i — Acos(2az + 4> + ir/2) as approximations to the exact 
functions derived in section EOl 
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Figure 9. The 4 fitting parameters for the simultaneous optimization for OT1 and OT2 for azimuth and elevation error. 
Units are in arcmin, except for the phase which is in rad. The lower two panels show the deformation corrected azimuth 
and elevation pointing errors p az and p e i , respectively. 
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Figure 10. The local azimuth and elevation deformation errors for OT1 and OT2. Units are in arcmin. The bottom panels 
are the difference between the two local deformations, showing a similar structure as in the Figsf3]and[S] 



elevation pointing error, respectively. 

We then further define the relative efficiency, rj re i, per antenna pair (i.e. per correlation) as the square root 
of the product of their primary Gaussian beams with an overlap error S re i (primary synthesized beam). For 
the absolute array efficiency, rj a bs, additionally an array synthesized beam product with an overlap error 5 a bs 
is multiplied with r\ re i and then averaged over the total number of array baselines ribase- This means that the 
efficiency is specifically calculated for a source with an extension comparable to the synthesized beam: 



Vrel = r=r- [ \J A(r, a Pl 0) • A(r, cr P ,S err ) d 3 r, (13) 

JdV 

^ n base 11/" I" 

Vabs = V — • — • / y/A{r,a P ,0) ■ A(r, crp, 8 e rr t )d 3 r ■ / A(r, a s , 0) • A(r, a s , S abSi ) d 3 r, (14) 

n baS e ^ N P N s J dv J dv 

where A(r, a, fx) = ^ ■ exp(— ^"^i ) is a Gaussian beam approximation, up and <r s are the primary antenna 
and synthesized beam FWHM, respectively. Np and N s are their normalization factors. 

In order to estimate the efficiency we run a Monte-Carlo simulation where we assume all the errors {d, m, p) 
to be uniformly distributed in their corresponding measured error interval^!. Fig llll illustrates the different 
efficiencies for an initial set of error intervals (in arcmin): d G [—1,1], m G [— 1, 1], p G [—1,1], each having an 
azimuth and an elevation component. For the following discussion we consider the error intervals for d and m to 
be fixed, d being a hard limit and m being considered an achievable limit after the radio alignment tuning. All 
the three errors can add up or compensate each other, where S a b s is more important when trying to improve the 



s Strictly speaking this assumption is not entirely correct because some of the errors are not independent, but correlated 
through the regular deformation patterns. 
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Figure 11. Efficiencies from a Monte-Carlo simulation for 10.000 realizations for an initial set of error intervals: d € [—1, 1], 
m £ [— 1, 1], p £ [—1,1]. Top panels: separated relative and absolute errors for one single baseline. Bottom panels: relative 
and absolute errors multiplied for a single baseline and for the entire array. The average overall efficiency is ~ 67%. The 
mean and median values are shown with a red and green marker, respectively. 

efficiency. S re i ~ 0.98 is not of a concern at all. In a scries of Monte-Carlo simulations ^ a fc s has been calculated 
for successively smaller error intervals for p, shown in Fig ll2l 



From the Figs fTTI and [T2l it is clear that the overall efficiency reaches a first limit, average and median ~ 0.67, 
when all the error components have the same error interval, which is set by d £ [—1,1]. However, successively 
smaller pointing errors can still improve the efficiency by about up to 10% due to the random summation and 
compensation of the errors. The final limit is reached at ~ 0.8 where smaller errors are absorbed and dominated 
by the largest error interval. A similar improvement (~10%) would be possible if m or d could be made smaller. 
An additional overall factor of about 0.9 will slightly reduce all these numbers, assuming that the distorted 
synthesized beam can be restored to about 90% with the phase correction model. 



4. SUMMARY AND CONCLUSION 

We have demonstrated that the local platform deformation error can be isolated from the pointing error with 
the help of two OTs. The extracted deformation error is in the range — V to +1' which is consistent with the 
photogrammetry data. The verification of an improved pointing accuracy is, however, non-trivial, because the 
platform deformation component will always be present on the ccd images. Therefore, one always has to rely 
on the modeling. Aiming for a pointing accuracy of ~ 0.5' rms (f] a b s ~ 0.76) seems reasonable. Even smaller 
pointing errors improve the efficiency only marginally. Similarly, the platform deformation induced phase error 
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Figure 12. Overall efficiencies for successively reduced pointing error intervals, from left to right, top to bottom: d 6 
[-1,4], [-1,3], [-0.5,0.5], [-0.2,0.2], [-0.1,0.1], [-0.05,0.05]. The mean and median values are shown with a red and 
green marker, respectively. 



is corrected with the photogrammetry data, and the synthesized beam amplitude is restored to 90% or more. 
Finally, one has to ask whether the separation of platform and pointing error is needed at all. Assuming that 
the errors are not separated and the derived pointing errors are then reduced to zero, one still ends up with 
an effective pointing error (p az and p e i in Eq. (fl"2|) ) in the range [—1,1] due to the over-/under-compensated 
platform deformation error. 6 a bs can therefore not be further reduced and the efficiency is limited to 67%. For 
large pointing errors (> 1') the separation becomes unnecessary. Going beyond rj a b s ~ 0.67 requires a better 
effective pointing and a separation of the error components. In any case, a phase correction model is needed to 
restore the synthesized beam. 
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